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Evidence for mixing in the neutral D meson system has recently been reported. Assuming negligible CP violation, 
non-vanishing width and mass differences between the two neutral D mass eigenstates has been found. Theoretical 
predictions of these are rather difficult, obscuring detection of New Physics contributions. However, the observation 
of CP violation in the D system would be a good signal of New Physics. We briefly describe the formalism that 
describes the neutral D decay and mixing, and present a method to determine all the mixing parameters accurately 
"qq ' allowing for arbitrary CP violation. 
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1. INTRODUCTION 

O . 

Neutral meson mixing is a flavor changing process, resulting in a flavor change AF = 2 and within the Standard 
£ — . 1 Model it occurs via box diagrams with internal quarks and W bosons. While the charged current flavor violating 
CI interactions can appear at the tree level, flavor changing neutral current (FCNC) interactions are feasible only at the 

, , loop level. Hence these interactions have an important role to play in the search for Physics beyond the Standard 

Model (SM), or New Physics(NP), as the new particles could appear virtually in the loops. One hopes to constrain 
I ■ NP by a measurement of these FCNC processes. 

Untill March 2007, neutral meson mixing had only been seen in the down type mesons: A"(in 1956), Bd(m 1987) 
and -B s (in 2006). The parameter €k in K° — K° mixing played a constraining role in building of models of NP. The 
measured mixing in the Bd and B s mesons, being consistent with SM predictions has also resulted in constraining 
t-H ■ various NP models. D meson is the only up type meson where mixing is possible. Evidence for mixing in the neutral 
D meson system has recently been reported by the Belle, BaBar and CDF collaborations. These 

experiments hnd non-vanishing width and mass differences between the two neutral D mass eigenstates assuming 
0\ ■ negligible CP violation (CPV). The HFAG average for ICHEP08 rules out the no mixing scenario at 9.8<t. In the 
Standard model one expects the mixing parameters in the D system to be small. Further, decays of the D meson 
via tree diagrams as well as D° — D° mixing- due to the negligible contribution of the internal b quark in the box 
diagram, essentially involve only two generations. Hence one expects that there should be no CPV in the charm 
■ system. 
> ' 

•i—l , 

2. THEORETICAL ESTIMATES AND LOOKING FOR NEW PHYSICS 

& ' 
_ _ i 

While the internal charm quark makes the dominant contribution in the box diagram for K° — K° mixing, the 
large top mass is responsible for appreciable mixing in Bd and B s mesons inspite of the suppression due to the small 
CKM elements. In fact, mixing in the neutral K and Bd systems resulted in predictions for the charm and top 
quark masses respectively, before direct discovery. However, in the D meson system, the heaviest down type quark, 
the b quark is not heavy enough to compensate for the large suppression due to the small CKM factor, IK^V^I 2 . 
Contributions from the internal light quarks (s or d) are dominant and hence, the mixing parameters in the D system 
are expected to be small. In the flavor SU(3) limit one would expect mixing to exactly vanish. In fact, it has been 
shown to arise only at second order in SU(3) breaking [7j. 

An explicit calculation of the SM mixing parameters in the D system is very hard. This is due to the fact that the 
D meson mass lies in the intermediate range where neither the inclusive nor the exclusive approaches work too well. 
An inclusive approach is based on an operator product expansion in terms of matrix elements of local operators of 
increasing dimensions with coefficients in powers of A/m c . However, m c is not large enough to allow the expansion 
upto few terms to be accurate. Such calculations [l0( suggest mass and width differences of order 10 -3 . On the 
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other hand, in an exclusive approach, summing over hadronic final states also fails as D is not light enough, that 
just a sum over few exclusive channels could suffice. Moreover, cancellations between states within a given SU(3) 
multiplet requires that the contribution of each state be known with high precision. In Ref. H SU(3), breaking 
arising from phase space differences was studied and it was concluded that the width difference could be at the level 
of one percent. 

Since the theoretical predictions of the D mixing parameters by various groups vary over orders of magnitude, a 
clear signal of NP from comparison of the expected and observed values of the mass and width differences may be 
hard. However, an observation of CPV will imply presence of NP [8j, independent of hadronic uncertainties. In most 
extensions of SM, the decay amplitudes and width difference are not expected to be affected by NP, however, the 
mass difference could be modified by new short distance CP violating contributions. In our discussion, we consider 
CPV only in the mixing and no direct violation Q . 



3. FORMALISM AND NOTATION 

Time evolution of the D system is governed by the Schrodinger equation: 
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where Mn = M 22 , r u = T 22 from CPT invariance. The neutral D mass eigenstates are related to the weak 

q\ 2 = 1 is obtained by diagonalizing 



eigenstates by, \D\ t ; 



q 
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\ Mia - i^f 

the Hamiltonian matrix. The off diagonal elements of the mass matrix are due to transitions via off-shell intermediate 
states while those of the decay matrix are from on shell states and therefore constitute respectively the the dispersive 
and absorptive parts of D mixing. If the magnitude of q/p differs from unity and/or the weak phase <f> — aig{q/p) is 
nonvanishing, this would signal CP violation. 

The time evolution of the states \D°(t)) and \D°(t)) which start of as pure \D°) and \D°) at t — is given by 
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and M and T are the average mass and width of the two mass eigenstates, also the the mass and width differences 
of these eigenstates are popularly written [Tl| in terms of the dimensionless variables, 
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r r 

The time evolution functions f±(t) are studied in the limit x <C 1, y <C 1 and tT <C 1. In this limit we have 
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Hence, the time dependent decay rates for a D° decaying to a final state / and D° — > / and D° — » / have the 
form: 



\A(D°(t) ^/)| 2 = e-" [X f + YfTt + Z f (Ttf +■ 
\A(D°(t) -/)| a = e- r * + ijr* + Z, (r;) 2 +• 



(8) 
(9) 



4. MEASURING THE MIXING PARAMETERS 

This time dependence is utilized in determining the mixing parameters. At t = 0, the only term in the amplitude 
of decay of D° is the direct amplitude D° — ► /. At any time t > 0, there is a mixing contribution through the 
sequence D° — > D — » /. The interference of this mixing contribution with the direct decay, involves the mixing 
parameters: x, y, 0, as well as the magnitude r and strong phase S of the ratio of the D° — > / and D° — ► / 

amplitudes and plays a key role in their measurement. 

The branching ratios of the Cabibbo favored (CF) decays are large and naively one might expect that one should 
use these decays to determine the parameters. However, if D° — > / is CF, D° — > / is Doubly Cabibbo Suppressed 
(DCS) and hence the interference term is too tiny compared to leading term in the time dependent decay amplitude. 
For example, the decay rate for D° K~tt + is, 



r(£>° - K-n+) = \A Kv \ 
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y cos(<5 - 4>)] Tt 



where, Ak-k = A(D° — > K ir + ) and the ratio of the DCS to CF amplitude is defined as: 
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Since rk-n , x and y are much less than unity, even the linear term in Tt is negligible compared to the constant term 
and the CF decay rate can can only be used to determine |Ak"tt|- In a DCS mode on the other hand the constant, 
linear and quadratic terms in Tt in the time dependent decay rate are all of the same order, allowing all three terms 
to be measurable. The decay rate for the DCS mode D° — > K + ir~ is, 
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where, due to the presence of a relative strong phase between the DCS and CF amplitudes, the combinations: 
x 'k* — (^coSfW+J/sintfifc) and y' Kv — (ycosS K „ — x sin S K7I ) appear. In the linear and quadratic terms, the suppression 
from the small mixing parameters is compensated by the larger D° — > f rate. The time dependent decay rate of 
D° —> K + ir~ as well as its conjugate mode was used by BaBar and CDF collaborations to determine [2j, |j| y' (from 
the linear term) and x' 2 (from the quadratic term), assuming CP conservation. 

In the case of singly Cabibbo suppressed (SCS) CP eigenstates modes, the strong phase is identically zero; and 
hence, the time dependent decay rate for these modes, like D° — > K~K + reduces to the simple form: 



T{D° 



K+K-) = \A KK \ 2 e~ ri 



1 - I - 1 (Vkk cos 4> - Xkk sin <j))Tt 
P 



Unlike the DCS modes where the term quadratic in Tt is enhanced by the ratio of CF to DCS rates, in the SCS modes 
all time dependent terms are of the same order in sin# c , hence quadratic and higher terms in Tt cannot be extracted. 
Assuming \q/p\ ~ 1 and = 0, the linear term in Tt can directly measure y, as has been done in Ref. [lj] (yep °1 
Ref. Q is y for no CPV). However, the time dependent study of only the SCS CP eigenstates does not allow x to be 
determined, even in the limit \q/p\ ~ 1 and = 0. 

A Dalitz plot analysis [3[ of D° — > K s ir + ir~ has also been performed by the Belle collaboration to determine all 
the mixing parameters. But this has systematic errors associated with the parameterization of the resonant content 
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of the Dalitz plot and hence is model-dependent. Measurements of yep using D° — > K S K + K~ as well as that of y' 
and x' 2 using D° —* K + tt~it were also reported [l2| at this conference. 

Since the SM estimates of x and y are uncertain, the current experimental results for these parameters cannot have 
any clear implications for NP. However, these results already constrain the parameter space of various models. In 
supersymmetric models with quark-squark alignment, constraints on the up- type squark matrices have been discussed 
in Ref. 13j. They seem to imply squark and gluino masses above 2 TeV. A detailed analysis of various NP models 
has been carried out in Ref. [141 ] . Out of the 21 models considered, only 4 are ineffective in producing charm mixing 
at the observed level. For the rest of the 17 models constraints on masses and mixing parameters are obtained. 

As pointed out earlier, within the SM, CPV in the D system is negligible and an observation of CPV would be a 
clear signal of New Physics. Babar and Belle have looked for CPV by calculating y' and x' 2 for D° and D° separately. 
They have also searched for CPV by measuring the difference of the decay rates of D° and D° to SCS CP eigenstates. 
No evidence of CPV has been obtained. It is hence important to have a technique to accurately measure the CP 
violating phase. 



5. DETERMINATION OF THE CP VIOLATING PHASE ALONG WITH OTHER MIXING 
PARAMETERS 

A technique to accurately determine all the mixing parameters including the CP violating phase has been given 
in Ref. 15|. It was shown that using the DCS mode D — ► K* a ir° and its conjugate modes, one can solve for all the 
D — D mixing parameters. This is possible if the K*° / K*° is reconstructed both in the self tagging if ± 7r =F mode 
and in the CP eigenstate K s tt° mode. 

With the K*° /K*° reconstructed in the self tagging K^ir^ modes, the time dependent decay rate has a form similar 
to that in Eq. (fT2"|) . The amplitude \A K ^\ = \A(D° — > K*°tt°)\ can easily be measured using the time integrated rate 
for the CF mode D° — > K* tt . The magnitude of the ratio of the DCS to CF amplitude can hence be determined 
using constant term in time dependent decay rate. The quadratic terms in Ft in the time dependent decay rates of 
D° — > K*°tt° and its conjugate mode D° — > K*°ir° readily determine \q/p\ and x' 2 + y' 2 . The interference terms, 
now involve the 4 unknown parameters: x, y, (j) and S, with a known value of x 12 + y' 2 . To determine them all, an 
additional observable is required. 

The K*° /K*° in the final state could also have been reconstructed in the neutral K s ir° mode, this provides the 
required additional observable. A unique feature of the final state K s ir°ir° is that it includes contributions from both 
K*°ir° as well as K*°ir° states; the amplitude for this final state is thus a sum of the CF and DCS amplitudes, 

|Ac a „| 2 = \A(D° -> K s tt\°)\ 2 = \A^\ 2 (l + r^ -2^ cos***) . (13) 

Since the decay mode involves two neutral pions it will not be easy to perform a time dependent study. Hence, we 
consider only the time integrated decay rate for this mode. The amplitudes A(D° — ► K s ir°ir°) and A(D° — ► i^ s 7r 7r°) 
are equal since i^ s 7r°7r° is a CP eigenstate. Hence, the time integrated decay rate for D° —> K s ir°ir° has the form, 

f \A(D°(t) -> K s TT°ir°)\ 2 dt w |yW| 2 [l +-(y cos<p~x sin <£) - 2r^ cos <W| . (14) 
Jo V 

Using this along with the linear terms of the time dependent decay rates of the self tagging modes allows a solution 
for tan 2 <f> and for x/y with a four- fold ambiguity, x and y can thus be individually determined since x' 2 + y' 2 is 
known. The solution obtained is finite even if cp — 0, with a correction term of order (cc cos 4rtr + ysbx8i(*„) sm<p. 
Hence an accurate estimation is possible, even if <j> is tiny. It should be possible to determine \x\, \y\ to order 7 x 10 -4 , 
4 x 10~ 4 respectively and <f> to about 1° at a Super-B factory with an integrated luminosity of 50 a6 -1 . 

While the SCS CP eigenstates like D — > K + K~ cannot alone be used to determine all the mixing parameters, 
minimal additional information from DCS modes makes this possible. This approach may provide the optimal method 
to determine all the parameters with current data. For the SCS-CP eigenstates, the strong phase is identically zero 
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and the ratio r = 1. The coefficients of the linear term in Tt, in the time dependent decay rate is a function of 3 
parameters: x, y and <fi. As pointed out earlier the quadratic and higher terms in Tt cannot be extracted. However, 
if we also include in this analysis the quadratic terms in Tt from the time dependent decay rates of DCS modes such 
as Kir, all the mixing parameters can be solved without approximation. Since x' 2 + y' 2 — x 2 + y 2 , these quadratic 
terms readily determine \q/p\ and f 2 = x 2 +y 2 . Alternatively, \q/p\ and f 2 could be measured using time integrated 
wrong sign relative to right sign SL decay rates. Having obtained \q/p\ and f 2 , 4> and x/y can easily be determined 
from D — > K + K~ and the solutions have been shown to be finite even for small 0. 

It may be further pointed out that if information from K + K~ modes is added to that from the K*~k modes it 
helps in reducing the ambiguities in x and y from four-fold to two-fold. It has recently been proposed to use the 
singly Cabibbo suppressed (SCS) D — > K*K modes to determine the mixing parameters 17]. However, if <j) is 
zero, these methods would be feasible only if the strong phase involved is measured elsewhere. 

6. CONCLUSIONS 

Mixing in the neutral D system has been clearly established. The Standard Model predictions involve large 
uncertainties, obscuring detection of New Physics contributions. An observation of CP violation in D mixing would 
clearly imply New Physics. The D — » K*°ir° modes are an example where it is possible to measure the CP violating 
parameters as well as the mass and width differences of the two D meson mass-eigenstates using only related final 
states, thereby reducing systematic errors. 
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